NASA-CR-200684 


/' v' - 


/' • ' -■ 

.J 

Ozone Production and Loss Rate Measurements in the Middle Stratosphere 
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The first simultaneous measurements of HO*, N0 r , and Cl x radicals in the middle strato- 
sphere show that NO* catalytic cycles dominate loss of ozone (O 3 ) for altitudes between 24 and 
38 km; Cl r catalytic cycles are measured to be less effective than previously expected; and there is 
no “ozone deficit” in the photochemically dominated altitude range from 31 and 38 km, contrary 
to some previous theoretical studies. 

Introduction 

Understanding the rate of removal of stratospheric ozone through catalytic cycles involving 
H0 X (HO 2 and OH), N0 X (NO 2 and NO), and halogen (CIO and BrO) radicals is essential for 
assessing the response of ozone to anthropogenic and natural perturbations such as industrial 
release of chlorofluorocarbons (CFCs) and halons, emission of nitrogen oxides from subsonic and 
supersonic aircraft, rising levels of N 2 0 and CH 4 , and enhanced levels of sulfate aerosols following 
volcanic eruptions. 1-3 In this paper we give the results of the first experiment to make simultaneous 
measurements of ozone-controlling radicals in the middle stratosphere, combine the results with an 
observationally constrained photochemical model, and show that production and loss of ozone are 
in photochemical equilibrium over much of this altitude region. 

Recent calculations from models of the mid-latitude stratosphere give the following theoretical 
results: catalytic cycles involving N0 X radicals are expected to dominate ozone destruction from 
about 22 to 45 km, while H0 X cycles dominate both above and below; 3-9 and modeled photochem- 
ical loss of ozone exceeds modeled production by 10 to 50% in the upper stratosphere. 4-8 The latter 
finding has been coined the “ozone deficit” problem, because at these altitudes ozone is expected 
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to be in photochemical equilibrium (production equalling loss), since the photochemical lifetime of 
ozone is much shorter than transport replacement times. 

The first finding has recently been tested below 21 km using in-situ aircraft measurements of 
HO x , N0 X and Cl x radicals . 10 These measurements, made under conditions of high aerosol loading 
as a result of the eruption of Mt. Pinatubo, are in good agreement with models which include the 
heterogeneous hydrolysis reaction of N 2 O 5 on sulfate aerosols . 3,4 

The ozone deficit problem above 45 km has also been recently discussed , 9 using satellite mea- 
surements of several long lived species, with the result that estimated production and loss of ozone 
are in better agreement than had been found in previous analyses of other satellite data. However, 
this study lacks direct measurements of H0 X or Cl x radicals. 

In this paper, we use the the first simultaneous measurement of concentrations of H0 X , N0 X , 
Cl x , and radical precursors throughout the middle stratosphere 11 to calculate the photochemical 
removal removal rate of ozone between 20 and 38 km. This analysis is the first to use direct and 
simultaneous measurements of radicals to determine the ordering of the loss cycles of ozone in the 
middle stratosphere, and to study the ozone deficit problem. 

Diurnally averaged photochemical removal rates of O 3 by each family of radicals are deter- 
mined in two ways: from direct measurements of radical concentrations (radical method), and 
from radical concentrations calculated by a photochemical model constrained by measurements of 
radical precursors (precursor method). The radical method is used to derive, for the first time, 
empirical measures of the total loss rate of ozone as well as the relative removal rate by each 
family of radicals which regulate the abundance of ozone in the middle stratosphere (24-40 km). 
The precursor method provides a basis for comparing the empirical loss rates to theoretical rates 
constrained by appropriate environmental conditions, represented by the observed concentrations 
of long lived radical precursors. Empirical loss rates found using the radical method demonstrate 
that catalytic cycles involving NO x radicals dominate photochemical removal of 0 3 in the middle 
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stratosphere, consistent with theoretical predictions using the precursor method. Our measure- 
ments demonstrate that production and loss rates of ozone (found using both methods) balance to 
within the uncertainty of measurement (10%) for altitudes between 31 and 38 km. This finding, 
which reconciles a long standing discrepancy in stratospheric chemistry, is due primarily to our 
observation of a lower abundance of CIO than calculated by models that allow for production of 
HC1 only from the reaction of Cl and hydrocarbons . 11 

Measurements and Model 

Descriptions of the measurement technique and analysis methods 12 have been published pre- 
viously; thus only a brief summary is given here. Data were acquired on a balloon flight launched 
from Ft. Sumner, New Mexico (34°N, 104°W) on 26-27 September, 1989. Concentrations of CIO 
were measured at mid-day between altitudes of 20 and 40 km with the Jet Propulsion Laboratory 
(JPL) Balloon Microwave Limb Sounder (BMLS), which senses microwave emission from the limb. 
Concentration profiles of OH, HO 2 , H 2 O 2 , HOC1, HC1, NO 2 , HNO 3 , N 2 O, H 2 O, O 3 , and temper- 
ature were measured throughout a diurnal cycle with the Smithsonian Astrophysical Observatory 
(SAO) far-infrared Fourier transform spectrometer (FIRS- 2 ), which detects atmospheric thermal 
emission in a limb viewing geometry. 

The constrained photochemical model used in our analysis balances, over a 24 hour period, the 
production and loss of 35 reactive gases for the temperature, pressure, and latitude at which the 
observations were obtained . 13 Standard reaction rates and absorption cross sections were used , 14 
except that photolysis cross sections for H 2 0 2 , HN0 2 , and HN0 4 were extrapolated to longer 
wavelengths to accurately represent photolysis of these molecules at large solar zenith angles. A 
reaction probability of 0.1 was used for the heterogeneous hydrolysis of N 2 O 5. 14 The heterogeneous 
hydrolysis of CINO 3 and BrNC >3 were included , 15 but have a negligible effect on model results for 
the temperatures of these observations. The altitude profile for the surface area density of sulfate 
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aerosols was adopted from SAGE II extinction measurements for this time period. 1 Other inputs to 
the constrained photochemical model include the profiles of temperature, O3, and H2O measured 
by FIRS-2. Profiles of CH 4 , odd nitrogen (N 0 y = NO + NO2 + NO3 + 2 xN 2 0 5 + HN 0 2 + 
HNO3 + HNO4 + CINO3 + BrNOs), inorganic chlorine (Cl y = HC 1 + C 1 N 0 3 + CIO + H 0 C 1 + 
Cl + 0 C 10 + 2XCI2O2 + C 100 ), and inorganic bromine (Br y = HBr + Br0N02 + BrO + HOBr 
+ BrCl) have been inferred from the FIRS-2 measurement of N 2 0 using relations derived from 
previous satellite 16 and in situ observations. 17 The maximum Br y at high altitude is assumed to 
be 21 ppt. 


O3 Photochemical Loss and Production Rates 

Ozone and atomic oxygen (O), together defined as odd oxygen ( 0 X ), are in rapid photo- 
chemical steady state via the reactions 03-fhi' — ^ 0+02 and 0 + 0 2 +M— ►Os+M. Throughout the 
stratospheric altitude range discussed in this paper, O3 represents more than 99% of total odd 
oxygen, thus a change in the abundance of odd oxygen is considered as equivalent to a change in 
the abundance of O3. Oxides of hydrogen, nitrogen, chlorine, and bromine regulate the abundance 
of 0 X in the middle stratosphere through a variety of reactions that catalyze the self reaction of 
O3 or the reaction of O3 and atomic oxygen. Since each reaction sequence listed below accounts 
for the loss of two molecules of 0 X , the net loss rates are determined by doubling the diurnally 
averaged rate of the limiting step (indicated by the number label) of each cycle. 

The role of H 0 X radicals in regulating the abundance of odd oxygen was first noted by Bates 
and Nicolet. 18 The important cycles are: 


OH + 0 3 -> H 0 2 + 0 2 , 

HO2 + O3 — ► OH + O2 + O2, (1) 

Net : 20 3 — 30 2 ; 
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and 

0H + 0 3 ^H0 2 + 0 2 , 

H0 2 + 0^0H + 0 2 , (2) 

Net : 0 3 + 0 -*■ 20 2 ; 

with other cycles accounting for less than 5% of ozone loss by H0 X species between 35 and 40 km, 
and less than 1% below 35 km. 

Crutzen 19 and Johnston 20 established the importance of removal of odd oxygen by cycles 
involving NO x : 

NO + O 3 — N 0 2 + 0 2 , 

N0 2 + 0 -4 NO + 0 2 , (3) 

Net : 0 3 + 0 — ► 20 2 ; 

and 

NO + 0 3 — + N0 2 + 0 2 , 

N0 2 + 0 3 -NO 3 +O 2 , (4) 

N0 3 hv — * NO + 0 2 , 

Net : 20 3 — * 30 2 - 

Molina and Rowland 21 and Wofsy et al. 22 drew attention to removal of odd oxygen by chlorine 
and bromine radicals, respectively. Loss occurs primarily through the cycles: 

Z0 + 0-+Z+0 2 , (5) 

z + o 3 — + zo + o 2 , 

Net : O 3 + 0 — * 20 2 ; 
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and 


H0 2 + Z0 H0Z + 0 2 , (6) 

HOZ + hv -+OK + Z, 

Z + O3 — ♦ zo + 0 2 , 

OH + Oa - H0 2 +0 2 , 

Net : 2O3 — * 30 2 ; 

where Z = Cl or Br, as well as through the cycle 

BrO + CIO -Br + Cl + Oz, (7) 

Br + O3 — * BrO + 0 2 , 

Cl + 0 3 -^C10 + 0 2 , 

Net : 20 3 -*• 30 2 ; 

which is a major loss cycle for O3 in the Antarctic spring. 

Loss of odd oxygen proceeds also by the reaction: 23 

0 + O3 — ► 20 2 ; (8) 

which is a minor sink in the middle stratosphere. Odd oxygen is produced by photolysis of 0 2 at 
wavelengths shorter than 200 nm, with a small contribution from the reaction of CH 3 0 2 and NO. 10 

The precursor method uses a photochemical model, constrained by our measurements of O3, 
H 2 0, N 2 0, and temperature, to predict the abundances of 0, H0 2 , OH, N0 2 , CIO, and BrO over 
a 24 hour period. These abundances are used to calculate the rates of reactions 1-8. The loss 
rates of 0 3 are integrated over 24 hours, using the DeMore et al. reaction rate constants, 14 and 
summed to evaluate the the rates of the HO x (1+2), NO x (3+4), halogen (5+6+7), and 0+0 3 (8) 
loss cycles. 
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The radical method uses direct measurements of HO 2 , OH, NO 2 , and CIO to calculate the 
rates of reactions 1-8. The measurements are interpolated onto a finer time grid by normalizing the 
abundances calculated by the constrained photochemical model to best fit the data. This process 
is shown in Figure 1 for a representative altitude of 32 km. We note that the model is in excellent 
agreement with the data even before scaling. We include 1 sigma uncertainties in the observations, 
plus any residuals between the observed and normalized concentrations. 

Figure 2 compares the removal rate of odd oxygen from the radical and precursor methods, 
partitioned into each radical family. This figure illustrates that the general behavior of the removal 
rates of odd oxygen inferred from the observed concentration of radicals is in broad agreement 
with theory for the H0 X and N0 X families of radicals. Theoretical values for the loss rate due to 
halogen cycles are calculated using two models: a conventional model which allows for production of 
HC1 only by reactions of Cl with hydrocarbons (a typical assumption of current multi-dimensional 
models of atmospheric photochemistry and transport 3 ’ 24 ’ 25 ); and a second model which includes 
additional production of HC1 from a 10% branch of the reaction of CIO with OH. DeMore et 
al . 14 suggests a lower limit of 0 %, and an upper limit of 14%, for production of HC1 from the 
CIO+OH reaction. Laboratory measurements of HC1 from this reaction, which would involve 
breaking and reformation of multiple chemical bonds, give ambiguous results. Our observations 
of CIO (see figure 1 ), supported by our simultaneous measurements of OH, HO 2 , HOC1, and 
HCl , 11 give best agreement with the the second model. Above 32 km, the lower abundance of 
CIO calculated by the second model gives removal rates of ozone by halogens which are reduced 
by as much as 40% compared to the conventional model. While this reaction is not the only 
possible pathway for additional production of HCl which might be consistent with the simultaneous 
measurement of CIO, HOC1 and HCl, it is the only pathway within suggested uncertainties of 
current kinetics measurements . 14 Corroborating evidence for the existence of such an additional 
pathway for production of HCl is provided by ATMOS measurements of HCl and CINO 3 6 , and 
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UARS measurements of HC1, CIO, and CINO 3. 26 

The fractional contribution of each radical family to loss of odd oxygen is shown in Figure 
3. Cycles involving N0 X radicals dominate the loss of ozone between altitudes of 24 to 38 km, 
contributing a maximum of 66 % to the total loss rate at 32 km. In contrast to observations 
obtained during SPADE , 10 our measurements also suggest that loss of odd oxygen due to N0 X 
cycles is comparable to loss due to halogen cycles near 20 km. The SPADE measurements showed 
loss of odd oxygen due to NO x cycles to be a factor of 2-4 slower than loss due to halogen and HO x 
cycles. However, the large uncertainties of our measurements at these altitudes make it difficult to 
determine precisely the ordering of the loss cycles. The balloon-borne observations were obtained 
in 1989, when concentrations of Cl y were approximately 14% lower than contemporary levels 27 and 
sulfate aerosol loading was about a factor of 5 less than values encountered during SPADE, which 
were elevated due to the eruption of Mt. Pinatubo in the summer of 1991. Both factors contribute 
to a larger relative contribution by NO x cycles to odd oxygen removal in in the lower stratosphere 
during 1989. 

Figure 4 illustrates diurnally averaged total loss rates for odd oxygen inferred from both the 
radical and precursor methods, as well as production rates from photolysis of O 2 calculated using 
a radiative transfer model that includes Rayleigh and aerosol scattering . 6,28 Three conclusions can 
be drawn from this figure. First* odd oxygen loss rates inferred from the observed concentrations 
of radicals (radical method) agree with theoretical loss rates (precursor method) to within the 
measurement uncertainty at all altitudes. Second, between 31 and 38 km, where the photochemical 
lifetime of ozone is short, we find that production and loss rates of odd oxygen calculated using 
both methods balance to within the uncertainty of the measurements of the radicals (approximately 
10 %), in contrast to previous studies (which lacked simultaneous measurements of HO 2 , NO 2 , and 
CIO) that reported loss rates up to 50% greater than production rates . 4-8 Third, rates of production 
below 31 km exceed the loss rates calculated using both methods by 30 to 40%, indicating that 
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this region is a net photochemical source region for O3. This imbalance is consistent with results 
of 2-dimensional photochemical-transport models that show production and subsequent transport 
of O3 to higher latitudes from this region. 8,29 

Our finding that production and measured loss of of odd oxygen balance for altitudes between 
31 and 38 km differs from the results of previous studies for two reasons: (i) we observe lower 
concentrations of CIO at these altitudes than are predicted by models which assume that production 
of HC1 occurs only by reaction of Cl with hydrocarbons, 11 as illustrated in Figure 1; and (ii) we use 
a formulation for photolysis of O2 that results in deeper penetration of ultraviolet radiation, 6 and 
better agreement with measured transmittances, than is found using the formulation recommended 
by WMO. 30 Including the additional source of HC1 lowers the total modeled loss rate of odd 
oxygen by about 17% near 38 km resulting in close agreement between production and loss of ozone 
found using the precursor method. Our findings are consistent with the study of Minschwaner et 
a/., 6 who postulated, based on constraints imposed by Atmospheric Trace Molecule Spectroscopy 
Experiment (ATMOS) measurements of O3, NO, NO2, N2O5, HNO3, CINO3, HC1, H 2 0, and CH 4 , 
that production and loss of odd oxygen would balance near 40 km if a minor branch of the reaction 
of CIO with OH results in production of HC1. Similar conclusions were reached by Chandra et 
a/., 31 who found that submillimeter heterodyne measurements of CIO and HC1 32 and the observed 
abundance and seasonal variation of O3 near 40 km were better simulated by allowing for the 
additional pathway for production of HC1. A more recent theoretical study 9 on the balance of 
production and loss of odd oxygen based on constraints provided by HALOE measurements of 
HC1, H2O, and O3 found a balance between production and loss between 35 and 40 km and excess 
production above 45 km. 
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Discussion 

Most current models which are used to assess the impact of anthropogenic emission of halogens 
on ozone overestimate the abundance of CIO in the upper stratosphere by neglecting the possibility 
of an additional pathway for production of HC1. 24 Our observations demonstrate that these mod- 
els overestimate the relative contribution of halogen cycles to the removal rate of odd oxygen at 
altitudes near 40 km. The predicted decrease in ozone near 40 km resulting from the build up of 
halogens during the past several decades is mitigated by 50% if an additional pathway for produc- 
tion of HC1 is included in models, resulting in closer agreement between observed and predicted 
trends in O3 at these altitudes. 31 Above 45 km, HO x catalyzed cycles dominate ozone loss, with 
the importance of halogen cycles decreasing with higher altitude into the mesosphere. Inclusion 
of an additional source of HC1 may not be sufficient to solve the imbalance between modeled pro- 
duction and loss at altitudes where reaction 2 dominates loss of odd oxygen. 5 The recent study of 
Crutzen et al. 9 based on HALOE measurements of HC1, H2O and O3 suggests excess production 
of ozone above 45 km. However, Crutzen et al. relied on theoretical profiles for OH and H0 2 ■ New 
observations of HO x radicals will be necessary to advance our understanding of the budget of O3 
at altitudes higher than 38 km. 

Since the loss rates found using the radical method depend on measured radical and O3 concen- 
trations and not modeled values for radical concentrations, the uncertainties in the calculated loss 
rates are dependent primarily on the measurement uncertainties and rate constant uncertainties in 
the rate limiting steps. At the altitudes where we measure a balance between production and loss 
of odd oxygen, 32 and 38 km, the dominant odd oxygen loss cycles involve reactions 2,3,5 and 8, 
with reaction 3 contributing more than 60% to the odd oxygen loss. Reaction 5 is dominated by the 
CIO reaction, which is more than 2 orders of magnitude more effective than the BrO reaction at de- 
stroying odd oxygen above 30 km. As a result, uncertainties caused by using modeled BrO at these 
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altitudes are very small with respect to other errors being considered. Using the average measured 
temperature of 240 K for these altitudes, the uncertainties in the rate constants for reactions 2,3,5, 
and 8 are 30%, 21%, 27% and 40%, respectively. By weighting these uncertainties by the relative 
contributions from each of these cycles, we calculate a net uncertainty in the calculated total loss 
rate of 23%. This value is on the order of the estimated uncertainty (20%) in the production of 
O3 6 and is larger than the uncertainties (10%) in the loss rates calculated from the measurements 
of the radicals. 

The accuracy of this set of radical measurements in the middle stratosphere, combined with 
the consistency between loss rates found using the radical and precursor methods, demonstrates 
that our first order understanding of processes that regulate the partitioning of radicals is correct, 
and that models, provided they allow for production of HC1 either by the reaction of CIO and 
OH, or some other currently unknown process which results in similar chlorine partitioning, should 
provide a realistic description of the relative rate of removal of odd oxygen in the lower and middle 
stratosphere by each of the major catalytic cycles. 3 
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Figure Captions 

Fig. 1. Individual concentration measurements (solid circles) are shown along with the modeled 
diurnal concentrations (solid curves) and the corresponding scaled concentrations (dotted curves) 
for HO 2 , NO2 and CIO at 32 km. The 1 sigma uncertainties are dominated by random measure- 
ment errors but include estimated systematic errors. The modeled CIO is also shown for the case 
of no production of HC1 from the C10+0H reaction (dashed curve). For this altitude, the average 
integration time for each measurement is 22, 42, and 50 minutes for H0 2 , N0 2 , and CIO, respec- 
tively; each integration comprises several shorter observations spread over average corresponding 
intervals of 2, 5.5, and 7 hours. 

Fig. 2. HO x (solid triangles), NO x (solid squares), and halogen (solid circles) catalyzed loss rates of 
O3 inferred from radical method and calculated using the precursor method (curves). Two halogen 
catalyzed precursor method loss calculations are illustrated: one assuming HC1 production only 
from Cl reacting with hydrocarbons and one with additional production of HC1 from a 10% branch 
of the OH -(- CIO reaction (see text). The O+O3 reaction curve can be considered to represent both 
the radical and precursor methods since 0 is inferred directly from measurements of O3, which is 
the precursor for 0. The symbols are offset in altitude for clarity. 

Fig. 3. The fractional total loss of odd oxygen from the HO x , NO x , and halogen catalytic cycles is 
shown as a function of altitude, as inferred from the radical method using measurements of H0 2 , 
NO 2 , CIO and O 3 (symbols), and from the precursor method (curves). 

Fig. 4. Loss rates of ozone calculated with the radical method (see text) from direct measurements 
of H0 2 , N0 2 , CIO, and 0 3 (solid circles), loss rates calculated using the precursor method (dotted 
curve), and production rates of O 3 (solid curve) calculated primarily from photodissociation of O 2 
(see text) and constrained by measurements of O 3 and temperature. 
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